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1. Introduction 



Debye screening of the color charge, expected at high temperature, is signaled by screening 
masses/lengths that can in principle be obtained from the gluon propagator [p. More specifically, 
chromoelectric (resp. chromomagnetic) screening will be related to the longitudinal (resp. trans- 
verse) gluon propagator computed at momenta with null temporal component (soft modes). In 
particular, we expect the real-space longitudinal propagator to fall off exponentially at long dis- 
tances, defining a (real) electric screening mass, which can be calculated perturbatively to leading 
order. Also, according to the 3d adjoint-Higgs picture for dimensional reduction, we expect the 
transverse propagator to show a confining behavior at finite temperature, in association with a non- 
trivial magnetic mass (see e.g. [Q]). We note that these propagators are gauge-dependent quantities, 
and the (perturbative) prediction that the propagator poles should be gauge-independent must be 
checked, by considering different gauges. 

Even though the nonzero-r behavior just described has been verified for various gauges and 
established at high temperatures down to around twice the critical temperature T c [Q, [|], it is not 
clear how a screening mass would develop around T c . At the same time, lattice studies of Landau- 
gauge gluon propagators at finite temperature in pure SU (2) and SU (3) theory have observed a 
sharp peak in the infrared value of the electric propagator around the deconfinement temperature, 
suggesting an alternative order parameter for the QCD phase transition [f|, g, g, 0, §]. (Of course, 
a relevant question is, then, whether this singularity survives the inclusion of dynamical quarks 
in the theory In the following, we investigate the critical behavior of electric and magnetic 
gluon propagators and try to characterize the screening masses around the transition temperature 
T c by performing large-lattice simulations in pure SU(2) gauge theory. In particular, we use the 
knowledge gained in the study of the zero-temperature gluon propagator (see [10] for a review) to 
identify systematic effects in the infrared limit and to define temperature-dependent masses for the 
region around and below T c . A more detailed analysis of our data will be presented elsewhere [p"l|]. 
(Preliminary results were reported in Jl2| , p"3|].) 



2. Results 



We have considered the pure SU(2) case, with a standard Wilson action. For our runs we 
employ a cold start, performing a projection on positive-Polyakov-loop configurations. Also, gauge 
fixing is done using stochastic overrelaxation and the gluon dressing functions are normalized to 
1 at 2 GeV. We take /3 values in the scaling region and lattice sizes ranging from N s = 48 to 
192 and from N t = 2 to 16 lattice points, respectively along the spatial and along the temporal 
directions. Our (improved) procedure for determining the physical temperature T is described 
in [OJ]. The momentum-space expressions for the transverse and longitudinal gluon propagators 
Dj (p) and D^{p) can be found e.g. in We first describe our investigation of the critical behavior 
for transverse and longitudinal gluon propagators and then discuss our proposal for computing 
screening masses around T c . 

Our data for the transverse (magnetic) propagator Dj{p) at the critical temperature are shown 
in Fig. |l]. We clearly see the strong infrared suppression of the propagator, as expected, with a 
turnover at around 400 MeV. Regarding systematic errors, there are considerable finite-physical- 
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Figure 1: Transverse gluon propagator at T c , for various choices of lattice size and j3. Values for the 
temperature, xN t , /3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given in 
the plot labels. The solid lines are fits, described at the end of this section. 



size effects, evidenced by the fact that the lattices with the smaller physical spatial size (the red, 
yellow, magenta and black curves) show different behavior from the remaining curves (green, blue 
and cyan). The nature of these effects is similar to what is observed at T = [10]. We mention that 



essentially the same features are seen for Dj{p) at other temperatures, above and below T c [O]. 

The longitudinal (electric) propagator Di(p) at T c is shown in Fig. ||. We immediately see 
severe systematic effects for the smaller values of N t . Let us note that our runs were initially 
planned under the assumption that a temporal extent N t = 4 might be sufficient to observe the 
infrared behavior of the propagators and our goal was, then, to increase N s significantly, to check for 
finite-size effects. As seen in Fig. ^, this assumption is not justified for the longitudinal propagator 
around the critical temperature, especially in the case of larger N s . Indeed, as N s is doubled from 
48 to 96 and then to 192, we see that the infrared value of Di{p) changes drastically, resulting in a 
qualitatively different curve at N s = 192, apparently with a turnover in momentum. (Note that, in 
this case, the real-space longitudinal propagator manifestly violates reflection positivity.) We took 
this as an indication that our choice of N t = 4 was not valid and therefore considered larger values 
of N t . We assume here that data points at N t = 16 are essentially free from systematic effects, since 



(as shown in [13]) the curves for temperatures around T c stabilize for N t > 8. As seen in the figure, 
we obtain in this way a different picture for the critical behavior of D^p). It is interesting to note 
(see Fig. ^ that the jV s effects at T c are significant for N t = 6 (with opposite sign with respect to the 
N t =4 case) and are still present for N t = 8 (and maybe also for N t = 16). This is also true slightly 
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Figure 2: Longitudinal gluon propagator at and around T c , for various choices of lattice size and j3. Values 
for the temperature, xNt, /3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are 
given in the plot labels. The solid lines are fits, described at the end of this section. 

below T c , but not immediately above T c [|l3|]. 

In summary, the transverse propagator Dj{p) shows significant finite-physical-size effects at 
T c , while the longitudinal propagator Dl(p) is subject to two sources of systematic errors for small 
N t : "pure" small-A^ effects (associated with discretization errors) and strong dependence on the 
spatial lattice size N s at fixed N t , when this value of N t is smaller than 16. The latter effect was 
observed only at T < T c , whereas the former is present in a wider range of temperatures around T c 
(see below). For all investigated values of the temperature, Di(p) seems to reach a plateau at small 
momentum p, while Dj{p) is infrared-suppressed, with a turnover in momentum roughly around 
350 MeV for all T ^ 0. 

Considering the infrared plateau in Di(p) — which we estimate here by Dz,(0) — as a function 
of temperature, the value observed at T = increases as the temperature is switched on, drops 
significantly for T >T C and then shows a steady decrease. In Fig. |3| we show data for Dl{0) for all 
our runs on the left-hand side, and for the region around T c on the right. We group together results 
from runs using the same value of N t , and indicate them by the label "DLfJJV/'. The data points 
indicated with "sym" correspond to symmetric lattices, i.e. to the zero-temperature case. Note that 
results for different N s 's at fixed N t may not fall on top of each other, which gives us an indication 
of the systematic errors discussed above. These are especially serious for N t = 4 around T c (red 
points). We see that, surprisingly, the maximum value of Dl(0) is not attained for T = T c — as 
might have appeared to be the case from the N t =4 lattices only — and it does not correspond to a 
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Figure 3: Infrared-plateau value for the longitudinal gluon propagator [estimated by Dl(0)] as a function of 
the temperature for the full range ofT/T c values (left) and for the region around T c (right). Data points from 
runs at the same value of N t are grouped together and indicated by the label "DL0_A^", where "sym" is used 
to indicate symmetric lattices (i.e. T = 0). 



flat curve from 0.5 T c to T c , as could be expected by looking only at these two temperatures. Rather, 
the maximum seems to lie at about 0.9 T c . Moreover, it clearly corresponds to a finite peak, which 
does not turn into a divergence as N s is increased at fixed N t . 

We have also looked at the real-space propagators. We find clear violation of reflection posi- 
tivity for the transverse propagator at all temperatures. For the longitudinal propagator, positivity 
violation is observed unequivocally only at zero temperature and for a few cases around the crit- 
ical region, in association with the severe systematic errors discussed above. For all other cases, 
there is no violation within errors. Also, we always observe an oscillatory behavior, indicative of 
a complex-mass pole. Typical curves for the longitudinal and transverse propagators in real space 
are shown (for T = 0.25T C ) in Fig. f| 

We now address the problem of characterizing the screening masses around T c . For all fits 
shown in Figs, [j] and ||] we used the same expression, namely a five-parameter fitting form of the 
Gribov-Stingl @, |J] type 1 

This form allows for two (complex-conjugate) poles, with masses m 2 = a ± ib, where m = 
triR + imj. The mass m thus depends only on a, b and not on the normalization C. The parameter 
77 should be 1 if the fitting form also describes the large-momenta region (from our infrared data 
we get r\ / 1). Recall that at high temperatures one usually defines the electric screening mass 
as the scale determining the exponential decrease of the real-space propagator at large distances, 
which is equivalent to D^O) -1 / 2 in the case of a real pole. We therefore expect to observe mj — > 



'Note that, for given values of a, b, d, 7], the global constant C is fixed by the renormalization condition, so that 
there are only four free parameters in ( |2.l| ). 
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Figure 4: Longitudinal (left) and transverse (right) gluon propagator in real space for T = 0.25T C . Values 
for xN t , j3, lattice spacing a and spatial lattice size L (both in fm, in parentheses) are given in the plot 
labels. Note that the solid lines are not fits. 



(i.e. b — > 0) for the longitudinal gluon propagator at high temperature. Note that, if the propagator 
has the above form, then the screening mass defined by D L (0)- 1 / 2 = ^{a 2 + b 2 )/C mixes the 
complex and imaginary masses niR and mj and depends on the (a priori arbitrary) normalization C. 

We generally find good fits to the Gribov-Stingl form (including the full range of momenta), 
with nonzero real and imaginary parts of the pole masses in all cases. For the transverse propagator 
Dj(p), the masses m# and mj are of comparable size (around 0.6 and 0.4 GeV respectively). The 
same holds for Di(p), but in this case the relative size of the imaginary mass seems to decrease 
with increasing temperature. A detailed discussion of the associated masses iur, mj is postponed to 
a forthcoming study [ 1 1], as we are presently considering variants of the above fitting form inspired 
by zero-temperature studies. Indeed, the use of a Gribov-Stingl form is motivated by the behavior 
of the gluon propagator at T = 0, where this type of expression has been shown to describe well 
the data in three space-time dimensions [[If]]. Recently, in [|l7|], various fitting forms of this type 
were used to describe large-lattice data for the 3d and 4d gluon propagator at T = 0. Noting that 
the 3d case may be considered as the T — > °° limit of the 4d case, we propose to interpolate the 3d 
and 4d zero-temperature forms to describe our finite- T data in 4d. 



3. Conclusions 



We study the longitudinal (electric) and transverse (magnetic) gluon propagators in momentum 
space, proposing the calculation of screening masses through an Ansatz from the zero-temperature 
case. Going from zero to nonzero temperature, we see that the electric propagator Di(p) is en- 
hanced, with an apparent plateau value in the infrared, while the magnetic propagator Dj{p) gets 
progressively more infrared-suppressed, with a clear turnover in momentum at all nonzero tem- 
peratures considered. Severe systematic effects are observed for the electric propagator around 
T c , suggesting that lattices of temporal extent N t > 8 are needed for this study. Once these errors 
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are removed, the data support a finite maximum (located at about 0.9 T c ) for the infrared value of 
D L (p). 
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